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BVI photometry of the very old open cluster Berkeley 17* 
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ABSTRACT 

We have obtained BVI CCD imaging of Berkeley 17, an anticentre open cluster that com- 
petes with NGC 6791 as the oldest known open cluster. Using the synthetic colour magnitude 
diagrams (CMD) technique with three sets of evolutionary tracks we have determined that its 
age is 8.5 - 9.0 Gyr, it distance modulus is (m — M)o= 12.2, with a reddening of E{B — V) = 
0.62 - 0.60. Differential reddening, if present, is at the 5% level. All these values have been 
obtained using models with metallicity about half of solar (Z=0.008 or Z=0.01 depending on 
the stellar evolution tracks), which allows us to reproduce the features of the cluster CMD bet- 
ter than other metallicities. Finally, from the analysis of a nearby comparison field we think 
to have intercepted a portion of the disrupting Canis Major dwarf galaxy. 

Key words: Hertzsprung-Russell (HR) diagram - open clusters and associations: general - 
open clusters and associations: individual: Berkeley 17 
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1 INTRODUCTION 

It is generally accepted that open clusters (DCs) are among the best 
tracers of the prope rties, the fo rmation and the evolution of the disc 
of our Galaxy fe.g.. lFrieil99l . This is particularly true for the old- 
est ones, since they could trace back the first moments of the disc 
lifetime and put crucial constraints on its formation mechanism(s) 
and early evolution. 

Berkeley 17 is among the oldest DCs of the Galaxy, if not 
the oldest one. It is located at 02000 = 05''20'"32', 02000 = 
+30°33'47", corresponding to the Galactic anticenter direction 
/ = 175.68, b = —3.68. This clust er had not received much atten- 
tion since its discovery bv iSetteducati & Weavej until 1994, 
when the almost contemporary works bvlKaluznv 1 ^3) (hereafter 
K94) IPhelps. lanes. & Montgomery! jl994h and Ijanes & Phelp^ 
definitely found it to be a very old object. The controversy 
on the actual age of Be 17 began immediately, since K94 stated that 
Be 17 is at least as old as NGC 679 1, whose age he assume d to be 9 
Gyr, while lPhelps et alj jl994l) and llanes & Phelps! iT994l) derived 
an age, based on the Morphological Age Indicator (MAI, see Sect. 
3) of about 12 Gyr, i.e., similar to that of the youngest globular 
clusters (albeit with a large error bar). A very large age has been re- 



* Based on observations made with the Italian Telescopio Nazionale 
GaUleo (TNG) operated on the island of La Palma by the Fundacion Galileo 
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iterated bv lPhelpslfl99^ (hereafter P97 : 121:^ Gyr), but not fully 
confirmed b y several other anal y 
9 ± 1 GYr: lsalariseTai]l2004 
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confirmed b y several other anal yses (e.g.. lCaiTaro et alj l998|.|l999|: 

10 Gyn TKmsbe^^^Chabovea 
I2OO6I : 10 ± 1 Gyr). 

The very large age found by P97 would imply that the age 
distributions of open and globular clusters overlap, and that the disc 
began to form stars during the last phases of halo formation. On 
the contrary, if Be 17 is conclusively demonstrated to be younger 
than the youngest globular clusters (GCs), the existence of a hiatus 
between the formation of the halo/thick disc and the present thin 
disc would be reinforced. The latter is nowadays the most popular 
scenario for the formatio n of our Galaxy and is supported by many 
observational evidences iFreeman & Bland-Hawthoml2002h . 

Less dispute appears to exist about the metallicity of Be 17, 
that has been e stimated to be slightly lower than solar by 
iFriel et alj tOOA with [Fe/H]=-0.33 (a = 0.12) dex based 
on low resolution spectroscopy, and by iFriel et all i2005h . with 
[Fe/H]=— 0.10 {a — 0.09) dex, based on high resolution spec- 
troscopy of three giant stars. A slightly sub solar or at most nearly 
solar abundance is also found by ph otometric analyses, both i n the 
optical and IR bands (see. e.g., P97: ICarraro et alll998lfl999l) . 

Even if its age is reduced to a value of 9- 10 Gyr, Be 17 is a fun- 
damental cluster because it defines the old tail of the OC population 
together with NGC 6791, from which it differs significantly. We are 
building a sample of homogeneously studied DCs well distributed 
in age, position and metallicity, th at we intend to use to study the 
disc properties and evolution; see iBragaglia & Tosil (2006) for a 
detailed description of our program (that we named Bologna Open 
Cluster Chemical Evolution - or BOCCE - program) and its results. 
Such a sample cannot leave Bel7 aside. We have obtained deep 
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Figure 1. Approximate positions of our pointings on Be 17 and the control 
field (larger squares); the smaller squares indicate the positions of the cen- 
tral field in K94 and P97. The map is 15 X 50 arcmin'^ , has North to the top 
and East to the left. 



B,V,I photometry that we present here and high resolution spec- 
tra of 3 red clump stars that will be the subject of a forthcoming 
paper. Medium resolution spectra of stars from the tip of the red 
giant branch to the main sequence are planned to determine mem- 
bership of more than 100 stars in this crucial cluster. 

The paper is organized as follows: we present our data in Sec- 
tion 2, we discuss the resulting CMDs in Section 3 and derive the 
fundamental cluster parameters using the synthetic CMDs method 
in Section 4; finally a discussion and a summary are given in Sec- 
tion 5. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations of Be 17 and of the associated control field used to 
check the field stars contamination were obtained at the Telesco- 
pio Nazionale Galileo on UT 2000 October 1, 2, and 3. We used 
DOLORES (Device Optimized for the LOw RESolution), a fo- 
cal reducer capable of imaging and low resolution spectroscopy, 
mounting a 2k Loral thinned and back-illuminated CCD, with scale 
of 0.275 arcsec/pix, and a field of view 9.4 x 9.4 arcmin^. The 
two fields, separated by about 30 arcmin, are shown in Fig.Q The 
observations include several exposures in each of the three BVI 
Johnson-Cousins' filters. Table 1 lists the dates of the observations 
together with the filters used and the corresponding ranges of expo- 
sure time (in seconds). All the three nights were photometric, with 
a mean seeing value around 1.4". 

Corrections for bias and flat-field were performed using the 
standard IRAF^ procedures. The subsequent data reduction and 
analysis was done following the same procedure for both data-sets 
(cluster and comparison field), and using the DAOPHOT II pack- 
age in IRAF (Stetson 1987; Davis 1994). Objects have been inde- 
pendently searched in all the frames, using a threshold of 4 ct above 
the background, and have been measured with a quadratically vary- 
ing point spread function (PSF). For each filter, the magnitudes 
have been aligned to that of a reference frame, the deepest one, 
obtained in the best seeing conditions. An average instrumental 
magnitude (weighted with the photometric errors) was derived in 
each filter. For the brightest stars, which are saturated in all frames 
but the shortest ones, only one measure was available. Finally, we 
applied aperture corrections (of the order of 0.2 - 0.3 mag) to the 
weighted instrumental PSF magnitudes in each band, to calibrate 
them on the same system of the standard stars. For the latter, mag- 
nitudes were derived using aperture photometry. 

The final catalogs have been created including all the objects 
identified in at least two filters, after applying a moderate selection 
in the shape-defining parameter sharpness (—2 ^ sharpness ^ 
2) and on the goodness-of-fit estimator (x^ ^ 10). To the two 
final catalogs, one for the cluster and one for the comparison field, 
we applied the transformation to astrometrize the a and S coor- 
dinates, using software written by P. Montegriffo at the Bologna 
Observatory. 



2.1 Photometric calibration 

Transformation from the instrumental to the standard system have 
been obtained using the standards areas PG023 1-1-051 and MarkA 
(Landolt 1992), observed several times during the three nights. 
We derived the calibration equations using only the standard fields 
observed on the same nights of the reference frames (October 
2 and 3). The 10 stars retained in our calibrations have colours 
-0.329 < B -V < 1.448 and -0.534 < V - I < 1.951, 
thus covering the bulk of the cluster stars without need for extrap- 
olation. 

For the extinction coefficients we used the average between 
the values of these two nights^ {kb ~ 0.23, kv' ~ 0.13, k/ = 
0.03). The resulting calibration equations have the following form; 



^ IRAF is distributed by the National Optical Astronomical Observatory, 
which are operated by the Association of Universities for Research in As- 
tronomy, under contract with the National Science Foundation 
^ see www.ast.cam.ac.uk/~dwe/SRF/camc_extinction.html and iKinJ 
ll985ft 
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Table 1. Log of observations for the cluster (02000 =05'' 20™ 32^52000 = 30" 33' 47"; / = 175.68, b = -3.68) and the external field (02000 =05'' 20' 
32**; (52000 = 31'' 03' 47"; I = 175.26, b = —3.40). Ny, Ng and N/ are the number of exposures (> 1) in each filter; t is the exposure time in seconds. 
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Figure 2. Comparison between our photometry and the one by K94 
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Figure 3. Comparison between our photometry and the one by P97 



B = b + 0.0460 X (& - v) + 1.4091 {rms = 0.013) 

V = v~- 0.1864 X (b-v) + 1.2310 (rms = 0.011) 

V = v- 0.0936 X {v-i) + 1.1895 (rms = 0.016) 

I = i + 0.0419 X {v~i)+ 0.7143 (rms = 0.013) 

where b, v, i, are the aperture corrected instrumental magnitudes, 
after correction also for extinction and exposure time, and B,V,I 
are the output magnitudes, calibrated to the Johnson-Cousins stan- 
dard system. 

Finally, we determined our completeness level in each band 
using extensive artificial stars experiments: we iteratively added, 
one at a time, about 50000 simulated stars to the deepest frames 
and repeated the reduction procedure, determining the ratio of re- 
covered over added stars (seelTosi et al. 2 0Q4. for a more detailed 
description). The results are given in Table 2. 

We checked the calibration comparing our photometry with 
that presented in previ ous literature (i.e., K94 and P97, which 
supersedes the one bv lPhelDS et"ai]ll994 taken in non optimal 
weather conditions). To estimate the differences, we cross identi- 
fied ou r objects with tho se of K94 and P97 (obtained through the 
BDA^. lMermillioJl995h . We found (see Fiesl2landl3l: (ABk) = 
-0.019, {AVk} = +0.030, {AIk) = 0.000 for K94, where A is 
intended as our magnitude minus the one by K94, and (ABp) — 
-0.121, (AVp) = -0.032, (Alp) = -0.024 for P97. Even if 
in neither case there are strong trends with magnitude, the compar- 
ison looks better with K94. The differences are small, except for 
our B photometry compared to P97's. The nights in which we ob- 
served the reference frames were deemed photometric, as was the 
case for the P97 observations (while K94 used the photometric part 
of a not completely photometric night), so we cannot easily explain 
this large difference. 

These differences in magnitudes translate into differences in 
colours: (A{B - V)k) = -0.049, {A{V - I)k) = +0.030 
and {A{B - V)p) = -0.089, (A{V - I)p) = -0.008: we are 
bluer than both literature photometries in the V, B — V plane; red- 
der than K94 and almost identical to P97 in the V,V — I plane. 
Since the simultaneous fit of the B — V and V — I colours is one 
of the indicators to select the most appropriate metallicity from the 
comparison with theoretical stellar models (see Section 4), such an 
uneven difference in the colours may affect the photometric metal- 
licity choice. 

We also compared the two photometries by K94 and P97 
with each other, finding differences in agreement with what is 
given above; if we limit the comparison to stars brighter than 

V = 20 as in P97 (their table 1), we find: {Bp - Bk) = 0.116, 
{Vp - Vk) = 0.060, and {Ip - Ik) = 0.058, based on about 
280 stars in common. This means {A{B — V)p-k) = 0.056, 



|http://www.univie.ac.at/webda/new.htmTl 
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{A( V - I) p-k) = 0.002. As already stated by P97, there are no 
direct means to decide which photometry is more precisely on the 
standard system (there are no stars with photoelectric measures in 
our fields), so we will proceed with our own data, without attempt- 
ing any correction. 

After all these comparisons were made , a new photometry has 
been p resented (but is not yet available) bv lKrusberg & Chaboveil 
yOOa). It appears of slightly lower quality than ours (e.g., the clus- 
ter sequences are less well defined) but, thanks to differential com- 
parison with NGC 188, it is claimed to be on the standard Landolt 
system to within ± 0.03 mag. Interestingly, they too find their pho- 
tometry bluer than the one by P97 by about 0. 1 mag in {B — V), 
providing further support to our calibration. 



3 THE COLOUR - MAGNITUDE DIAGRAM 

The final, calibrated sample of the cluster stars (that will be made 
available through the BDA) consists of 2473 objects identified in at 
least two filters, of which 1940 identified simultaneously in all the 
three filters BVI. 

The corresponding CMDs are shown in Fig.|4] where the mean 
errors per magnitude bin are also plotted. Both CMDs show: (i) a 
very clear main sequence (MS) extending down to V 24 for the 
V, F - / CMD and to V ~ 23 for the V,B-V CMD ; (ii) a main se- 
quence turn-off (MSTO) near V = 18, with a sparsely populated red 
giant branch (RGB), extending more than three magnitudes above 
the TO. This morfology is typical of very old stellar systems for 
which the RGB raises very steeply and is very extended vertically, 
while the colour separation between the MSTO and the base of the 
RGB is smaller than in intermediate age clusters. The red clump 
(RC), including ~ 8 stars, is visible at V ~ 15 and V — / ~ 1.7, 
B — V ~ 1.6. This clump coincides with tha t identified by K94 and 
P97 an d, in our opinion, is also apparent in lKrusberg & Chaboved 
i2006h CMD, in spite of these authors' conclusion. 

We have measured the difference in magnitude between the 
RC and the MSTO: SV ~ 2.9 ± 0.1. This, using the so called MAI 
jjanes&Pheip3ll994 . implies an approximate age of 14.5 Gyr, 
but even the conservative error of 0. 1 mag in the 5V value means 
that the cluster could be more than 16 Gyr old, or as "young" as 
about 12 Gyr Since it is now generally accepted that this calibra- 
tion leads to ages systematically larger than those derived e.g, from 
isochrone fitting, we take this only as an indication of the very old 
age of the cluster We intend to re-calibrate the age-5V^ relation 
once the entire sample of OCs in our BOCCE program is examined 
and homogeneously dated through stellar evolutionary models. 

In both CMDs the sequence of binary stars is visible, although 
its identification is complicated by the presence of a substantial 
field star contamination. For this reason, we did not attempt to 
quantify the fraction of binary syst ems as we did for le ss con- 
taminated cluste rs (e.g., NGC 6253: iBragaglia et'Zlll99l Be 21: 
iTosi et alll99l) . 

Finally, we note that information on membership based o n 
radial velocity is available only for a few bright stars. .Scott et al,^ 
il995h observed about 30 stars (all brighter than V ~ 16) deter- 
mining velocities with a precision of about ± 10 km s~^ and their 
data were used by iFriel et all i2002h for the determination of the 
cluster metallicities. They found that only about half of them could 
be considered good candidate cluster members. They measured a 
mean velocity of —84 (cr = 11) km s"'^ and a mean metallicity 
[Fe/H 1=— 0.33 (a = 0.1 2) dex. Four of these stars were later stud- 
ied bv lpriel et all J2005h . who find that only three are actual cluster 



members and derive an average velocity of — 73.7 (cr = 0.8) km s 
and a mean [Fe/H]=-0.10 (cr = 0.09) dex. 



3.1 The control field 

To estimate the field stars contamination we used the data of the 
control field, which is far enough from the cluster to be safely con- 
sidered free of cluster members. Fig.|5]shows the resulting V,V — I 
and V,B — V CMDs for the objects identified in this field. The 
corresponding catalogue contains 1310 objects identified in at least 
two filters, of which 1103 are identified simultaneously in all the 
three filters. Comparison of Figs.|4]and|5|helps to understand which 
are the true cluster features. We also plot (see Fig.|6j radial CMDs, 
both for Be 17 (where the centre was defined using histograms of 
the star coordinates along the x and y axes) and the comparison 
field. The cluster is quite loose, but its features clearly stand out 
with respect to the field population in the central part. 

As in many other cases (see e.g., iD'Orazi et alj Eoo3 for 
Be 32), the CMD of the external field appears to be composed by 
(at least) two populations. As already noted by K94, there is a broad 
main sequence that crosses diagonally the diagram and can be ex- 
plained by the young population of the Perseus spiral arm. As a 
de monstration, in Fig.|7t b) we have plotted the Z=0.02 isochrones 
bv lBertelH et"ai]<1994li for ages of 0.80, 1.25, 2, and 3 Gyr, assum- 
ing the values usually attributed to the arm E{B — V) =0.5 and 
(m — Af)o=l 1.75 (i.e., a distance from the Sun of about 2.2 kpc). 

There is a second sequence, fainter and bluer than the cluster 
MS (noticeable also in the CMD of the field centered on Be 17). In 
Fig.0b) we have overplotted on this second sequence the Z=0.008 
and Z=0.004 isochrones with older ages (4, 5, 6.3, 9.9 Gyr) and as- 
suming larger distance and reddening. We find a reasonably good fit 
with - V) ~ 0.65, (m-A/)o~ 13.35, i.e., a distance of about 
4.7 kpc, if Z=0.004 and E{B - V) ~ 0.60, (m - M)o- 13.30, 
i.e., a distance of about 4.8 kpc, if Z=0.008. With these values, 
we cannot have intercepted the Monoceros (or Anticentre) Ring, 
a feature of the Galactic disc visible towards the extern al parts of 
our Galaxy (e.g., iNewberg et all200lllbata et all2003h but with a 
much larger dist ance from the Sun, of abo ut 10 kpc. This Ring h as 
been associated iMartin et alj2004l but see lMomanv et alj2004l for 
a different view) to an over density in star counts, interpreted as 
the remnant of the dwa rf galaxy Canis Major (CMa). Interestingly, 
iBellazzini et alj i2006h present the N-body model of the disruption 
of CMa, which shows a clear signature towards the anticentre, and 
at a distance of about 5 kpc from the Sun (see their fig. 11). We 
are possibly seeing this feature, so we have compared the CMD of 
CMa to the one of our control field. We took the data presented by 
IBellazzini et alj i2004h for CMa, taking into account the different 
areas, as well as the distance modulus and reddening given in that 
paper, i.e. E{B - V) = 0.10 and (m - M)o= 14.5. The resuU is 
shown in Fig0c); the two CMDs are comparable, and we also see 
the blue plume (for V magnitudes between about 17.5 and 19, and 
colours bluer than about B — V = 1). These comparisons suggest 
that we are actually seeing a small portion of the disrupting dwarf. 

Fi nallv. FiglTTd) sho ws the output of the Besanjon Galactic 
model iRobin et alj|2003|) in the direction of our control field, af- 
ter correction for our completeness function; the bright parts of the 
CMDs are similar, but they are a poor match fainter than about 
V = 20. This shows that we cannot explain our observed features 
in terms of simple Galactic field, but we need to take the anoma- 
lous overdensities into account to reproduce all the sequences of 
the CMD. 
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Figure 4. Left panel: V,V — I CMD for Be 17; Right panel: V,B — V CMD. The mean errors per interval of magnitude V are also plotted. 
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Figure 5. Left panel: V,V — I CMD for the external field; Right panel: V,B — V CMD. The mean errors per interval of magnitude V are also plotted. 
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Figure 6. Radial CMDs for Be 17 (upper panels) and corresponding areas of the comparison field (lower panels); the center was assumed at pixel x=1225, 
y=l 175 in both cases, and we plot stars within distances of 1, 2, 3 arcmin from it. The CMDs contains 103, 385, 790 objects for panels (a), (b), (c) respectively, 
and 28, 177, 372 for panels (d), (e), (f) respectively. 



4 CLUSTER PARAMETERS 

Age, distance and reddening of Be 17 have been derived with the 
sam e procedure applied to a ll the clusters of our BOCCE project 
(see lBragaglia&Tosill200d and references th erein), namely th e 
synthetic CMD method originally described bv lXosi et alj il99lh . 
The best values of the parameters are found by selecting the cases 
providing synthetic CMDs with morphology, colours, number of 
stars in the various evolutionary phases and luminosity functions 
(LPs) in better agreement with the observational ones. 

As usual, to test the effects of the adopted stellar evolution 
models on the derived parameters, we have run the simulations 
with three different types of stellar tracks, with various prescrip- 
tions for the treatment of convection and of overshooting from con- 
vective regions. Actually, Be 17 is so old that its MSTO stars have 
masses around 1 Mq and are therefore not expected to experience 
any overshooting from convective cores. Indeed, we find no appre- 
ciable difference in the results obtained with models with or with- 
out overshooting. 

To estimate the metallicity which better reproduces the pho- 



tometric properties of the cluster, we have created the syn- 
thetic CMDs adopting, for each type of stellar models, metal- 
licities ra nging from solar down to 20% of solar. We still as- 
sume fsee lBragaglia & Tosil2006l) as solar metallicity models those 
with Z=0.02, because they are the ones calibrated by their au- 
thors on the Sun, independently of the circumstance that nowa- 
days the actual solar metallicity is supposed to be lower (see 
lAsplund. Grevesse. & Sauvaill2005h . At any rate, we consider the 
metallicities obtained with our photometric studies only indicative 
and use, whenever available, high resolution spectroscopy for a safe 
determination of the chemical abundances. 

The adopted sets of stellar tracks are listed in Table 3, where 
the corresponding references are also given, as well as the model 
metallicity and the information on their corresponding overshoot- 
ing assumptions. The transformations from the theoretical lumi- 
nosity and effective temperature to the Johnson-Cousins magni- 
tudes and colours have been performed using Bessel, Castelli & 
Pletz (1998) conversion tables and assuming EiV — I) = 1.25 
E{B - V) (Dean et al. 1978) for all sets of models. Hence, the 
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Figure 7. What kind of pupulation is present in the comparison field (panel a) ? In panel (b) we show that it can be the combination of a young, nearby 
component, with E{B — V) ~0.5, (m — J\f)o~11.75), corresponding to the Perseus arm, and an older, farther, metal-poorer one. The reddening and distance 
[E{B — V) ~0. 60-0. 65, (m — A/)o~13.30-13.35] of the latter suggest that we could have intercepted the orbit of the disrupting CMa. The younger population 
is fitted with the four brighter dark isochrones (red in the colour electronic version), with age 0.8, 1.25, 2 and 3 Gyr. The older population is compared with 
the four lighter (green in the colour version) Z=0.004 isochrones , and wi th the four dark (blue) Z=0.008 ones, corresponding in both cases to ages of 4, 5, 6.3 
and 9.9 Gyr. In panel (c) we plot the CMa CMP i Bellazzini et a ll2004) assuming the reddening and distance given above [after correction for t he values of 
E(B- V)=O.\0, fm-M)n=14.5 in lBellazzini et alji2004il and scaling for the different areas. Panel (d) shows the result of the Galactic model bv lRobin et alj 
i2003l) in the direction of our field, after correction for the completeness function. 



Table 2. Completeness level for the central and external fields; mag is the 
calibrated magnitude (B,V or I, calibrated with the equations given in the 
text, and assuming b — v and v — i = I). 



Table 3. Stellar evolution models adopted for the synthetic CMDs. The FST 
models actually adopted here are an updated version of the published ones 
(Ventura, private communication) 



mag 


cs 




cj 


cs 


cv 


cj 


16.00 


1.0 


1.0 


1.0 


1.00 


1.00 


1.00 


16.50 


1.0 


1.0 


0.97 


1.00 


1.00 


1.00 


17.00 


1.0 


1.0 


0.99 


1.00 


1.00 


1.00 


17.50 


1.0 


1.0 


0.96 


1.00 


1.00 


1.00 


18.00 


1.0 


0.97 


0.98 


1.00 


1.00 


0.98 


18.50 


1.0 


0.98 


0.96 


1.00 


1.00 


0.99 


19.00 


0.99 


0.99 


0.94 


0.98 


0.98 


0.96 


19.50 


0.99 


0.97 


0.88 


0.97 


1.00 


0.96 


20.00 


0.98 


0.97 


0.81 


0.98 


0.97 


0.92 


20.50 


0.98 


0.97 


0.63 


0.96 


0.97 


0.80 


21.00 


0.97 


0.95 


0.31 


0.96 


0.95 


0.48 


21.50 


0.94 


0.93 


0.13 


0.92 


0.97 


0.12 


22.00 


0.94 


0.91 


0.03 


0.89 


0.93 


0.02 


22.50 


0.91 


0.73 


0.0 


0.87 


0.87 


0.0 


23.00 


0.81 


0.45 


0.0 


0.66 


0.55 


0.0 


23.50 


0.57 


0.18 


0.0 


0.38 


0.29 


0.0 


24.00 


0.24 


0.04 


0.0 


0.13 


0.04 


0.0 


24.50 


0.05 


0.0 


0.0 


0.02 


0.0 


0.0 


25.00 


0.0 


0.0 


0.0 


0.00 


0.0 


0.0 



Set 


metallicity 


overshooting 


Reference 


BBC 


0.02 


yes 


Bressan et al. 1993 


BBC 


0.008 


yes 


Fagotto et al. 1994 


BBC 


0.004 


yes 


Fagotto et al. 1994 


FRA 


0.02 


no 


Dominguez et al. 1999 


FRA 


0.01 


no 


Dominguez et al. 1999 


FRA 


0.006 


no 


Dominguez et al. 1999 


FST 


0.02 


?j=0.00 


Ventura et al. 1998 


FST 


0.02 


r?=0.02 


Ventura et al. 1998 


FST 


0.02 


r?=0.03 


Ventura et al. 1998 


FST 


0.01 


r?=0.00 


Ventura et al. 1998 


FST 


0.01 


r?=0.02 


Ventura et al. 1998 


FST 


0.01 


?7=0.03 


Ventura et al. 1998 


FST 


0.006 


?7=0.00 


Ventura et al. 1998 


FST 


0.006 


?7=0.02 


Ventura et al. 1998 


FST 


0.006 


r?=0.03 


Ventura et al. 1998 



different results obtained with different stellar models must be as- 
cribed fully to the models themselves and not to the photometric 
conversions. 

The synthetic stars are attributed the photometric error derived 
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from the artificial stars tests performed on the actual images. They 
are retained in (or excluded from) the synthetic CMD according to 
the photometry completeness factors listed in Table 2. All the syn- 
thetic CMDs have been computed either assuming that all the clus- 
ter stars are single objects or that a fraction of t hem are members 
of bin ary systems with random mass ratio (see iBragaglia & Tosil 
I2OO6I for a description of how binaries are included in the synthetic 
CMDs). We find, as in many other clusters, that a binary fraction 
around 30% well reproduces the observed distribution along the 
main sequence. All the synthetic CMDs shown in the figures as- 
sume this fraction of binaries. 

In spite of the field contamination, the evolutionary sequences 
are quite well defined in the CMD of Be 17. Fig.|^shows that the 
CMDs of the inner cluster regions are not much cleaner than the 
global one, and have the disadvantage of containing less stars; we 
have thus run the simulations for the whole field covered by our im- 
ages. Since the cluster field contains 1940 stars measured 'mB,V,I 
and the external field 1 103, we assume that the cluster members are 
1940 - 1103 = 837. Hence, the synthetic CMDs have been created 
with 837 objects. 

We find that in all cases a solar metallicity must be excluded, 
because it does not allow us to simultaneously reproduce both the 
observed B ~ V and V — I colours. With Z=0.02, when the syn- 
thetic B — is correct, V — I always turns out bluer than observed. 
Vice versa, models with Z<0.006 provide V — I colours redder 
than observed, when B — V is. correct. On the other hand, stellar 
tracks with metallicity about half of solar can reproduce very well 
the observed colours of all the evolutionary phases, once the appro- 
priate reddening is adopted. For the BBC models, the only avail- 
able metallicity leading to self-consistent results is Z=0.008. For 
both the FRA and the FST models, Z=0.006 allows to reproduce 
rather well the observed colours, but Z=0.01 leads to a much better 
agreement. We thus favour Z=0.01 as the photometric metallicity 
of Be 17. 

Fig.[8|compares three representative cases of synthetic CMDs 
(panels b, c and d), with the empirical one (panel a) in the V, B — V 
plane. To better appreciate the differences between the various the- 
oretical predictions, in panels b, c and d we plot only the 837 
synthetic stars attributed to Be 17. Fig. |9| shows the correspond- 
ing V,V — I plane, with the superposition of the 1103 external 
field stars to the 837 synthetic ones. Hence, all the CMDs of Fig.|9| 
contain 1940 objects. The observed V ~ I cluster sequence looks 
slightly wider than the synthetic one, perhaps due to the character 
of the field contamination, but we think this does not significantly 
affect the age determination. It is also immediately apparent that 
our synthetic MSs never reach as deep as the empirical one. This 
is due to the fact that none of the adopted sets of stellar evolution 
models includes stars less massive than 0.6 M0, while our photom- 
etry clearly detects significantly smaller stars. 

The binary sequence is evident in the synthetic CMDs of Fig.|8| 
and it is interesting to notice that it perfectly overlaps with existing 
sequences of the observational diagram, one on the right of the MS 
and one above the MSTO. These two observational sequences are 
not statistically significant by themselves, due to the field contami- 
nation, but the perfect match with the predicted binary sequence is 
unlikely to be just a fortuitous coincidence. 

With all the adopted sets of stellar evolution models, the age 
that allows us to better reproduce all the observed properties of 
Be 17's CMD is between 8.5 Gyr (with the BBC models, panel b 
in Figs[8|and|9) and 9 Gyr (with the FST models, panel d in Figs 
[8|and|9). Namely, with this age we obtain the right luminosity and 
colours of both the MSTO and the clump, the right shape of the 
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Figure 8. Comparison between empirical and synthetic CMDs. Panel a 
shows the B — y CMD of the 1940 stars with measured B, V and I. Pan- 
els b, c and d show representative cases of synthetic diagrams of the 837 
assumed cluster stai's. Panel b illustrates the best case for the BBC models: 
it assumes Z=0.008, age=8.5 Gyr, E{B - V) =0.62, and (m - J\/)o=12.2. 
Panel c shows the BBC predictions for an age of 12 Gyr: the best metal- 
licity is still 0.008, E{B - V) becomes 0.56 and (m - M)o= 12.1. 
Panel d illustrates the best case for the FST models: Z=0.01, age=9 Gyr, 
E{B - V) =0.60 and (m - M)o=12.2. 
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Figure 9. Compaiison between empirical and syntlietic CMDs. Panel a 
shows the V — I CMD of the 1940 stars with measured B, V and I. Panels 
b, c and d show the same representative cases of synthetic diagrams of Fig. 
Isl but with the addition in the CMD of the 1 103 stars measured in the exter- 
nal field, to allow for a more direct comparison with the empirical diagram. 



MS, the subgiant branch (SGB) and the RGB, the right number of 
stars predicted on the clump and on the RGB, and the right star 
counts at the various MS levels. Younger ages imply fainter clump, 
worse MSTO morphology and, often, more extended SGB, and are 
therefore easy to reject. On the other hand, it is admittedly diffi- 
cult to completely rule out older ages, up to 11-12 Gyr. As shown 
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Figure 10. Luminosity functions of the 1940 stars in the cluster region as 
measured from the data (dots) and from the synthetic CMDs of Fig.|9] Error 
bars on the empirical values are within the size of the symbols. The solid 
line coiTesponds to the objects in panel b, the dotted line to those in panel c 
and the dashed line to those in panel d: they are barely distinguishable from 
each other. 



in panel c of Figs [H and |9l (BBC models with Z=0.008, age = 12 
Gyr, E{B - V) = 0.56 and (m - M)o= 12.1), the clump does 
become brighter and bluer, the star distribution on the RGB does 
differ somehow from the observed one, the SGB luminosity does 
increase too much with decreasing temperature, but none of these 
defects can really be considered bad enough to reject the case, once 
both the theoretical and the observational uncertainties are taken 
into account. However, since all the adopted sets of stellar evolu- 
tion models agree in favouring an age of 9.0 ±0.5 Gyr, we consider 
this value as the most likely age of Be 17. 

The various sets of stellar models converge to values in strik- 
ingly good agreement also in the predictions for the reddening and 
the distance modulus. We obtain E{B - V) = 0.62, (m - M)o= 
12.2 for the BBC models with Z=0.008 and age = 8.5 Gyr; E{B - 
V) = 0.60, (m - M)o= 12.2 for the FRA models with Z=0.01 and 
age = 9.0 Gyr; E{B - V) = 0.60, (m - M)o= 12.2 for the FST 
models with Z=0.01 and age = 9.0 Gyr. For sake of completeness 
(and because it was suggested by K94) we have also tested the pos- 
sibility of some differential reddening affecting this rather contam- 
inated region. From the comparison of synthetic CMDs computed 
with varying amounts of differential reddening, we suggest that, 
if present, the absorption variations should be rather small, with 
AE{B ±0.03. 

The LFs of the three synthetic CMDs shown in Fig.|9|are plot- 
ted in Fig. llOK lines) and compared to that of the CMD of Fig.|9^. 
The latter is significantly affected by the 1103 external stars, thus 
making the comparison of little help in discriminating between dif- 
ferent cases. We simply find that all these models are consistent 
with the data, except for the lack of very low mass stars which lead 
us to under predict the faint end of the LF. 



5 SUMMARY AND DISCUSSION 
5.1 Comparison with previous results 

As already mentioned in the Introduction, Be 17 has received a lot 
of attention in the last decade. 

The first to present good CMDs for this cluster was K94; he 
observed the central part of Be 17 with the 2.1 m telescope and a 
larger area (~ 23 arcmin^) with the 0.9 m telescope, both at Kitt 
Peak. He deduced a cluster diameter of about 13 arcmin and men- 
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tioned the possibility of differential reddening. By comparison of 
the CMD of the central part to NGC 6791 and assuming that Be 17 
is metal-poorer, he found an age equal or slightly larger than that 
of NGC 6791, an apparent distance modulus (?n — M)v > 15.0, 
and a reddening E{V - I) > 0.7 (or E{B - V) > 0.56). He also 
concluded that the field stars appearing redder and brighter than 
the cluster MS belong to the Perseus spiral arm. The results for 
age and reddening are in very good agreement with our findings. 
We find instead a much shorter distance, but cannot readily explain 
such a difference. 

The same year Ijanes & Phelpj il994) . in their study of 
the propertie s of the old OCs, using photometry presented by 
IPhelps et alj il994;) and their calibration of the MAI, derived the 
very large age of 12 Gyr for B e 17. The MAI had not b een intended 
to measure absolute ages, and ljanes & PhelpsHl994 stress that it 
is to be used only to measure relative ages. In their calibration of 
the 5V versus age (the latter taken from internally precise, but in- 
homogene ous sources) they also c onsidered several GCs, with ages 
taken from lChabover et al]il992h . Among them the youngest GCs 
have ages of about 10 Gyr. The MAI was scaled so that an age of 
15 Gyr was obtained for the oldest among the globulars: an age 
currently considered older than the Universe ! 

Another kind of analysis was done by P97, who concentrated 
on the single cluster, instead of considering the bulk properties of 
the open cluster sample. He obtained deep and precise B,V,I pho- 
tometry of Be 17, using the 2.1 m Kitt Peak telescope on a photo- 
metric night. P97 noted, as we already said in Sect. 2. 1 , that his and 
K94 photometries do not perfectly agree, but he thought that this 
did not compromise his major conclusions. He used a control field 
to statistically subtract the field star contamination and employed 
the cleaned CMDs for his analysis. First he determined the cluster 
age by means of the SV and the SI (this is a difference in colour, 
that can be calibrated to derive a second estimate of SV): he found 
5V = 2.7 (instead of 2.8 as in lPhelps et 11111994 and from the 
MAI calibration obtained an age of 10.9 Gyr Later he derived the 
cluster age from isochrone fitting. Using the ^VandenBera (1985) 
isochrones (that do not reach the RC) he derived a be st fit age of 
12 Gyr with Z=0.01, while for the lBertelli et all Jl994 isochrones 
he cited a range of 10 to 12 Gyr, with a preference for Z=0.02. 
On the whole, he assumed an age of 12+2 Ciyr and a metallicity 
between solar and half of solar The reddening and distance he de- 
rived are perfectly consistent with ours: 0.52 ^ E{B — V)^ 0.68 
and (m — M)o = 12.15 ± 0.10. With these values, and adding 
the largish diameter and mass and the rather high radial velocity, 
P97 seems to think that Be 17 may be a transitional cluster, in be- 
tween open and globular ones; his main conclusion is that this clus- 
ter poses important constraints on the age of the disc and of the halo 
too, which is thought to have formed before. 

W hy did he obtain such a large age using the iBertelli et alj 
il994) isochrones, coming from the same set of models from which 
we obtain instead a best fit age of 8.5 Gyr and Z=0.008 ? With his 
data we would have made a choice different from his. In fact, judg- 
ing from his figures 7 and 8, the luminosity of the RC is actually 
better fit by the 10 Gyr is ochrone. Furthermor e, we would have 
given lower weight to the IVandenBerd Jl985h isochrones, since 
they miss the very important clue of fitting both the MSTO and 
the RC luminosities and colours. 

A co nfirmation of our choice comes from the work by 
ICarraro et a l. ( 1998J: they re-analy sed the data by P97 (and K94), 
using the new lGirardi et alj i200(j) isochrones (in preparation, at 
the time) and found that the possible combinations of parameters 
are: E{B - V) = 0.55 - 0.67, (m - M)o = 12.3 ± 0.07, and 



age= 9 ± 1 Gyr for Z= 0.007 - 0.013, with the metallicity not 
affecting the a ge. They too d iffer from P97's choice of the best fit 
from the Bertell i et al. NI994 isochrone, and place lower weight to 
thelVandenBerg < 19851) ones. 

Another check of the cluster parameters come from the work 
bv lCarraro et alj il99St) based on J, K IR data. They obtained an 
independent estimate of the cluster metallicity from the RGB slope, 
finding [Fe/H]~ —0.35, in very good agree ment with the spec- 
troscopic estimates. Using this value and the lGirardi et alj i2OO0l) 
isochrones they find a best fit age of 9 Gyr; translating from the 
IR quantities, they derive E{B - V) ^ 0.58 and (m - M)o = 
12.0 ± 0.1. 

ISalaris etaljj2004 derived a calibration of the SV that takes 
into account also the metallicity. They determined distances to 10 
clusters by means of main sequence fitting, and the corresponding 
ages using their own evolutionary models. These clusters were used 
to derive the calibration from which, using literature values for SV 
and metallicity, an age of 10 Gyr is obtained for Be 17. They also 
confirm a gap of 2.0±1.5 Gyr between the formation of the thin 
disc and of the halo, while the oldest OCs (Be 17 and NGC 6791) 
have the sam e age of the oldest thin disc st ars. 

Finally, Krusberg & Chaboved i200d) assumed the metallic- 
ity from Friel et al. ( 2 002) and obtained, from isochrone fit with- 
out RC, an age of 10 ± 1 Gyr, (m - M)v=14.1-14.2 and a 
E[B — V) =0.56-0.61. The y do mention tha t assuming a higher 
metallicity (as indicated bv iFriel et alj|200'^ the age would de- 
crease, thus presumably becoming closer to our value. 

Since we have not analysed yet the spectra we obtained for 
three RC stars, no real comparison between our (approximate) pho- 
tometric metallicity based only on the best fit isoc hro nes and the 
spectr oscopic values obtained bv lFriel et al] j2002h and lFriel et all 
( 2005) is attempted. We only note that there is a general agreement 
that the cluster metallicity is (slightly) sub solar. 



5.2 Discussion 

Our analysis confirms that Be 17 is a very old disc object but 
that its ag e is far from tha t of globulars. There is only one known 
exception: ICarrard i2005h has presented evidence that GCs may 
be much younger than usually taken for granted. Whiting 1, with 
b — —60.6°, low metallicity (Z = 0.001), and heliocentric distance 
of 45 kpc hardly qualifies as a disc open cluster, yet it has an age 
of about 5 Gyr. Its nature has to be further investigated; either the 
separation between Galactic globular and open clusters is fuzzier 
than we think, or Whiting 1 is a "freak", maybe connected to an 
accretion event. 

Considering the age distribution of the bulk of well behaved 
globular clusters, recent studies of a ver y ample, well defined , 
homogeneously analy sed sample of GCs iRosenberg et aljfl999l : 
iDe Angeli etal]l2005l) have derived accurate relative ages for 55 
GCs, finding that most of them are coeval and old. There are some 
GCs younger by about 1 Gyr, but the really younger ones (by about 
3 Gyr) seem all connected to streams, i.e. not truly genuine halo or 
disc cl usters. The absolute a ge is more delicate to determine; for in- 
stance. lGratton et al]i2003h . after deriving a very accurate distance 
using the main sequence fitting method for three GCs also present 
in that sample, found that two of them, NGC 6752 and NGC 6397, 
are coeval while the metal- richer 47 Tuc i s youn ger by about 2.5 
Gyr. This, using models bv lStraniero et alj ^993), and taking into 
account differences in microscopic diffusion and helium content 
between those models and later results, means absolute ages of 
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about 13 and less than 1 1 Gyr, respectively, with an error bar of 
about 1 Gyr. 

Our result places Be 17 quite safely away from that lower 
limit; however, to obtain an absolute age we have to rely on the- 
oretical models, and ages of GCs and OCs obtained by different 
authors with different techniques and models may not be immedi- 
ately comparable. To get a definitive answer, stricter homogeneity 
is required, but the evidences are still largely in favour of the thin 
disc being younger than the halo. 

Finally, we note that a hiatus between the formation of the halo 
and thick disc and that of the presently observed thi n disc has also 
been f ound using chemical signatures. For instance. iGratton et alj 
i200(J) analysed disc and halo stars in the solar neighborhood and 
noticed that the run of [Fe/O] versus [O/H] supports the hypoth- 
esis of a period wi th no, or ve r y low, star formation, confirming 
previous results by iFurhrmam] il998h . In a [Fe/O] versus [O/H] 
plot (where O is the clock), Fe - primarily produced by type la Su- 
pernovae - after remaining constant with increasing O, suddenly 
increases at constant O. Since O is instead produced by short-lived, 
massive stars, this is interpreted as absence of such stars, i.e., of 
star formation, for a period of at least 1 Gyr. 

5.3 Summary 

We have analysed B, V, I photometric data of Be 17, and deter- 
mined its fundamental parameters by means of comparison of the 
observed CMDs to synthetic ones generated using different sets of 
stellar evolutionary models. 

• We find an age of 8.5-9 Gyr, a distance (m — M)o= 12.2 
(consistently for all evolutionary tracks used; an error of ± 0. 1 mag 
can be assumed), a reddening value E{B — V) = 0.60-0.62, and 
approximate metallicity Z=0.008 or 0.01. The three adopted sets of 
models agree very well on these values. 

• We cannot completely rule out a larger age ( ^ 12 Gyr) but 
the corresponding synthetic CMDs produce consistently worse fits 
to the observed ones. 

• Differential reddening is estimated to be at most at the 5 per 
cent level. 

• We find that a binary fraction of 30 per cent well reproduces 
the MS: these binary systems populate a secondary brighter and 
redder sequence that is possible to detect also in the observed 
CMDs along the single-stars MS and brighter than the MSTO. 

• The SV measured on our data would imply a very old age 
( > 10 Gyr), but a re calibration of this useful relative age indicator, 
based on homogeneously determined ages, is necessary. 

• Comparison of our findings with literature values indicates 
only two real discrepancies: with K94 for distance and with P97 
for age; possible causes for the latter have been discussed. 

• The comparison field shows a structure probably associated to 
the young population of the Perseus arm; there is a second compo- 
nent, compatible with having intercepted the orbit of the disrupting 
CMa galaxy. 

• Be 17 has a sub solar metallicity, a quite large reddening, and 
a very old age. With our preferred value of 8.5 - 9 Gyr Be 17 is 
definitely younger than the bulk of GCs, those that are thought to 
have formed in the very first phases of Galaxy formation, either by 
a monolithic collapse or early fragment accretion. 

Further work on Be 17 is foreseen, both to derive its elemen- 
tal abundanc es on a scale homog eneous with the other OCs in our 
sample (e.g.. lCarretta et all2005l and references therein) and to de- 
fine membership from radial velocities for about 100 stars in crucial 



evolutionary phases. The latter information may be used to confirm 
or refine the determination of cluster properties done in the present 
paper. 
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